DISTRIBUTION  STATEMENT  A:  Approved  for  public  release;  distribution  is  unlimited. 


Ultrashort  Laser  Pulse  Propagation  in  Water 

George  W.  Kattawar  and  Alexei  V.  Sokolov 
Dept,  of  Physics 
Texas  A&M  University 
College  Station,  TX  77843-4242 

phone:  (979)  845-1 180  fax:  (979)  845-2590  email:  kattawar@tamu.edu 
phone:  (979)  845-7733  fax:  (979)  845-2590  email  sokol@physics.tamu.edu 

Award  Number:  N000140810037 
http://people.physics.tamu.edu/trouble/work.html 


LONG-TERM  GOALS 

The  major  objectives  of  this  proposal  are  twofold.  The  first  objective  is  to  perfonn  both  an 
experimental  and  theoretical  study  of  the  factors  affecting  the  propagation  distance  behavior  of 
ultrashort  (femtosecond)  laser  pulses  in  water.  This  study  will  be  conducted  in  the  so-called  linear 
regime  which  involves  laser  intensities  that  are  below  the  threshold  where  nonlinear  effects  set  in.  A 
fundamental  problem,  which  will  be  resolved  by  this  research,  is  whether  or  not  the  temporal  width 
and  spacing  of  short  pulses  affects  its  absorption  spectrum.  If  one  simply  knows  the  spectral 
composition  of  the  input  pulse  and  the  absorption  spectrum  of  the  water,  is  this  sufficient  to  predict  the 
temporal  evolution  of  the  pulse?  If  such  is  not  the  case,  then  we  will  have  to  consider  the  way  the 
system  responds  to  pulse  widths  and  pulse  spacing  which  are  short  compared  to  vibrational  relaxation 
times  in  water. 

The  second  phase  of  the  research  will  be  to  explore  the  nonlinear  regime  where  dramatic  changes  to 
the  temporal,  spatial,  and  spectral  properties  of  the  medium  occur.  The  primary  processes  being  self- 
focusing  and  self-phase  modulation  due  to  the  Kerr  effect  (also  called  the  quadratic  electro-optic  effect 
which  was  discovered  by  John  Kerr  in  1875).  Self-focusing  can  lead  to  an  enormous  increase  in  the 
peak  intensity  where  long  filaments  can  occur  and  in  some  cases  lead  to  supercontinuum  generation  or 
“white  light”  generation  first  discovered  in  1970. 

OBJECTIVES 

The  propagation  of  an  ultrashort  pulse  of  light  through  a  linear  and  absorptive  medium  such  as  water, 
is  of  great  fundamental  importance  for  several  reasons.  One  of  the  most  important  of  which  is  that  it 
may  be  possible  to  transmit  information  over  much  greater  distances  using  ultrashort  pulses  compared 
to  propagation  distances  achieved  by  using  pulses  with  long  time  durations,  including  CW  (continuous 
waves).  The  first  prediction  of  optical  precursors  was  given  by  both  Sommerfeld  and  Brillouin1  in 
1914  using  an  asymptotic  method  now  called  steepest  descent.  Their  analysis  was  based  on  a  step- 
modulated  field  propagating  through  a  Lorentz  dielectric  which  is  nothing  more  than  a  collection  of 
damped  harmonic  oscillators.  Later  refinements  to  their  conclusions  were  made  by  Oughstun  and 
Sherman2.  The  first  measurements  which  claimed  to  observe  optical  precursors  in  deionized  water 
were  made  by  Choi  and  Osterberg3  where  they  found  that  the  precursors  were  attenuated  non- 
exponentially  with  distance.  What  was  extremely  significant  about  this  work  was  that  they  found  the 
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pulse  energy  detected  at  a  distance  of  5  m  in  water  with  a  60  fs  pulse  was  one  hundred  times  greater 
than  the  signal  obtained  using  a  900  fs  pulse  (having  the  same  fundamental  wavelength)  which  had  the 
usual  Lambert-Beer  behavior  ( It  should  be  noted  that  the  so  called  Lambert-Beer  law  was  first 
discovered  by  Pierre  Bouguer  in  his  published  work  in  1729  titled  Essai  d’optique  sur  la  gradation  de 
la  lumiere  and  we  will  hereafter  refer  to  it  as  the  Bouguer-Lambert-Beer  Law  or  BLB.  The  law  is  valid 
for  monochromatic  sources  but  unfortunately  researchers  have  mistakenly  applied  it  to  pulse 
propagation  where  the  bandwidth  covers  changes  in  the  absorption  coefficient ).  Their  conclusions 
were  brought  into  serious  doubt  by  Alfano4  et  al.  who  claimed  that  their  bandwidth  was  not  wide 
enough  to  cross  the  760  nm  absorption  band  of  water  and  therefore  that  their  conclusions  were 
questionable..  In  a  later  paper,  Fox  and  Osterberg5  found  deviation  from  BLB  when  they  used  pulses 
with  temporal  widths  of  60  fs  and  a  repetition  rate  of  1  kHz  where  they  observed  more  than  two  orders 
of  magnitude  less  absorption  after  propagation  through  6  m  of  water  when  compared  to  BLB.  They 
also  found  that  BLB  was  not  violated  for  pulses  with  varying  bandwidth  with  temporal  widths  of  900 
fs  (pulses  centered  at  800  nm)  and  repetition  rates  of  80  MHz.  Li  et  al6  have  performed  a  more  recent 
measurement  for  liquid  water  and  they  found  deviations  from  the  BLB  law  after  the  pulse  had 
propagated  only  1.5  m  and  the  deviation  continued  to  increase  out  to  3.5  m  which  was  as  far  as  they 
could  detect  the  signal.  The  shortest  pulse  they  used  was  10  fs  with  a  repetition  rate  of  75  MHz  which 
gives  a  pulse  spacing  of  13.3  ns.  They  showed  that  deviation  from  BLB  behavior  only  occurred  when 
the  pulse  bandwidth  was  of  the  order  of  100  nm  but  for  shorter  bandwidths  (<  30  nm)  BLB  behavior 
was  verified.  It  is  important  to  note  that  this  pulse  spacing  is  much  longer  than  the  vibrational 
relaxation  times  in  liquid  water  which  are  typically  of  the  order  of  a  picosecond.  The  most  recent 
experiments  by  Okawachi  et  al.7  used  four  different  laser  systems  centered  at  800  and  1530  nm  to 
show  the  dependence  of  pulse  duration  and  repetition  rate  on  pulse  absorption.  Their  bandwidths 
ranged  from  20-60  nm  and  they  showed  that  BLB  behavior  was  obeyed  in  every  case.  It  is  now  clear 
that  there  is  a  real  enigma  associated  with  the  entire  phenomena  of  ultrashort  laser  propagation  in 
water. 

APPROACH 

Experimental: 

There  are  three  basic  improvements  over  what  has  been  done  in  the  past  that  we  are  proposing: 

(a)  Improved  input  pulse  shapes 

The  laser  system  that  we  will  use  in  this  experiment  is  the  one  we  currently  use  for  femtosecond  CARS 
spectroscopy  of  dipicolinic  acid  and  bacterial  endospores8'11.  We  have  an  amplified  femtosecond  laser 
(Mira  +  Legend,  Coherent)  whereby  we  obtain  up  to  1  mJ  per  35  fs  pulse  (at  1  kHz  repetition  rate)  at 
the  fundamental  800  nm  wavelength.  We  can  control  the  chirp  (positive  or  negative)  of  these  800  nm 
pulses  over  a  wide  range  (increasing  the  pulse  duration  by  up  to  two  orders  of  magnitude,  to  over  3  ps, 
while  keeping  the  spectrum  unchanged)  very  easily,  and  very  precisely  by  adjusting  our  pulse 
compressor.  Thus  we  will  be  able  to  study  propagation  of  laser  pulses  of  different  length  but  having 
identical  spectra.  Fig.  1  shows  the  experimentally  measured  pulse  shapes  obtained  by  this  method  and 
measured  by  our  home-made  variation  of  frequency-resolved  optical  gating  (FROG). 
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Fig.  1.  Chirped  pulses  obtained  by  adjusting  the  pulse  compressor  and  measured  by  a  modified 
frequency-resolved  optical  gating  technique.  The  numbers  in  the  upper  right  corner  of  each  plot 
give  the  pulse  duration;  the  spectra  of  all  four  pulses  are  identical. 

Alternatively,  we  can  modify  the  pulse  spectral  intensities  while  keeping  the  spectral  phases  fixed. 

Our  home-made  pulse  shaper  allows  us  to  make  the  spectral  width  up  to  100  times  narrower  and  the 
pulse  length  correspondingly  longer  without  introducing  any  chirp.  Moreover,  our  pulse  shaper  allows 
us  to  fine-tune  the  precise  frequency  of  these  narrow-band  pulses.  We  will  use  such  spectrally- 
modified  pulses  for  the  detailed  study  of  light  propagation  through  water.  In  addition,  we  have  the 
capability  to  produce  some  more  complex  pulse  shapes,  should  it  become  necessary  or  interesting. 

Our  recent  paper12  gives  a  colorful  illustration  of  potential  possibilities,  and  shows  an  example  of  the 
use  of  pulse  shaping  for  coherent  Raman  spectroscopy.  More  complex  pulse  shapes  will  be 
particularly  important  for  the  studies  of  nonlinear  pulse  propagation.  In  the  second  year  of  the  project, 
we  plan  to  purchase  and  use  a  commercially  available  computer-programmable  pulse  shaper  for  fast 
(and  potentially  adaptive)  shaping  of  the  laser  pulses  propagating  through  water. 

Further  into  the  project,  we  will  expand  our  work  into  other  wavelength  regions,  and  use  our  two 
computer-controlled  optical  parametric  amplifiers  (OPerA,  Coherent),  pumped  by  our  amplified 
femtosecond  laser  system  (described  above).  Signal  and  idler  pulses  obtained  from  the  two  OP  As  can 
be  frequency-doubled  or  mixed  with  the  fundamental  to  produce  up  to  20  pj  per  50  fs  pulse  at  tunable 
visible  wavelengths.  In  addition,  we  may  use  un-amplified  light  from  the  fs  oscillator  (at  85  MHz),  for 
fast-repetition-rate  low-peak-power  experiments. 

While  the  previous  experimental  work  on  optical  precursors  suffered  from  limited  choice  of  laser 
sources,  our  laser  system  will  allow  great  flexibility  and  precise  control  over  the  pulse  shapes,  for  a 
broad  range  of  well-controlled  measurements. 

(b)  Spectrally-resolved  highly-sensitive  measurement  at  the  output 

We  will  measure  the  transmitted  spectral  intensities  (and  not  just  the  total  output  power).  We  will  use 
a  state-of-the-art  CCD-equipped  spectrometer  which  will  count  photons  at  each  wavelength  separately. 
A  second  CCD-equipped  near-infrared  spectrometer  will  allow  us  to  extend  our  overall  working 
wavelength  range.  What  has  been  used  in  many  of  the  past  experiments  were  photomultiplier  tubes 
(PMT’s).  A  PMT  is  a  "bucket"  detector  that  just  counts  the  total  number  of  photons  and  does  not 
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provide  spectral  resolution  (unless  it  is  combined  with  a  spectrometer,  in  which  case  a  tedious 
scanning  is  needed).  PMT’s  are  virtually  unsurpassed  in  sensitivity;  however,  a  good  cooled  CCD  can 
be  almost  as  sensitive.  Highly  sensitive  signal  detection  over  a  large  dynamic  range  is  essential  for  the 
success  of  the  proposed  project,  since  the  intensity  of  light  transmitted  through  many  meters  of  water 
will  vary  over  many  orders  of  magnitude  as  a  function  of  propagation  length.  Our  group  has  a  great 
deal  of  experience  in  performing  measurements  where  the  signal  varies  up  to  over  nine  decades.  As  an 
example,  we  reproduce  here  a  figure 13  [Fig.  2]  which  shows  signal  varying  over  seven  orders  of 
magnitude.  That  particular  set  of  data  was  obtained  without  spectral  resolution.  Currently  we  are 
capable  of  achieving  similar  or  even  higher  sensitivity  and  dynamic  range  with  spectral  resolution, 
which  is  afforded  by  the  use  of  a  spectrometer  equipped  with  a  cooled  CCD.  Our  measurements  on 
pulse  propagation  through  water  will  be  done  with  spectral  resolution. 


Fig.  2.  An  example  from  our  previous  work,  where  we  used  three  ultrashort  pulses 
at  different  wavelength  (810  nm,  660  nm,  and  575  nm)  and  measured  femtosecond 
coherent  anti-Stokes  Raman  scattering  (CARS)  signal  (measured  in  methanol-water  solutions) 

varying  in  magnitude  over  many  decades15 . 

At  a  further  stage  of  the  project  we  plan  to  measure  the  transmitted  pulse  shapes.  For  relatively  small 
propagation  distances  (several  meters)  we  will  use  a  commercial  autocorrelator  (presently  available  in 
our  lab),  and  a  FROG  apparatus  (that  we  intend  to  purchase).  For  large  propagation  distances  this 
measurement  will  become  challenging  since  the  output  power  will  get  very  low.  For  the  low-power 
measurements,  we  will  set  up  a  cross-correlation  measurement,  such  as  the  one  that  has  produced 
pictures  shown  in  Fig.  1.  Another  good  option  for  direct  pulse  shape  measurement  is  an  ultrafast 
streak  camera.  We  have  recently  installed  and  tested  such  a  state-of-the-art  streak-camera 
(Hamamatsu)  in  our  lab,  and  we  will  use  it  to  measure  the  duration  of  pulses  propagated  through  over 
10  meters  of  water. 

(c)  Water  cell  design 

We  have  designed  and  already  constructed  our  own  water  cell.  Our  focus  in  designing  the  new  cell 
was  to  minimize  the  number  of  mirror  reflections,  and  to  be  able  to  increase  or  decrease  the  path 
lengths  continuously  and  not  just  in  steps  of  0.5  meters. 
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In  previous  experiments,  the  maximum  path  length  through  the  water  was  6  meters.  In  our  water  cell, 
we  will  be  able  to  run  experiments  at  a  maximum  path  length  of  18  meters.  A  schematic  of  a 
simplified  experimental  set  up  is  provided  in  Fig.  3.  We  currently  have  the  equipment  installed  (as 
shown)  and  can  have  the  complete  setup  within  a  short  period  of  time. 


Specti  ornetef 


Fig.  3.  Experimental  setup  schematics. 


Theoretical: 

In  order  to  be  able  to  model  the  complete  pulse  propagation  process,  we  must  be  able  to  handle  the 
actual  pulse  shape  that  will  be  used  in  our  experiments  along  with  the  best  data  we  have  on  the  actual 
absorption  spectra  of  pure  water.  In  fact  the  water  we  will  use  will  be  prepared  by  the  same  method  as 
in  the  absorption  measurements  done  by  Pope  and  Fry14.  Since  we  have  a  great  deal  of  experience 
with  the  finite  difference  time  domain  (FDTD)  method,  we  will  have  to  substantially  rewrite  parts  of  it 
to  handle  the  actual  dispersion  of  the  refractive  index  of  water,  particularly  around  the  overtones  of  the 
OH-stretch/scissor  coupled  modes.  The  regions  around  the  sharp  slope  changes  in  refractive  index  are 
precisely  the  regions  where  the  interesting  precursors  will  appear.  This  type  of  approach  has  never 
been  performed  but  is  the  sine  qua  non  for  successful  modeling.  The  basis  of  our  method  is  to  take  the 
two  Maxwell  curl  equations  in  Rationalized  MKSA  units  and  assuming  no  currents;  namely, 
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and  then  discretize  them  over  a  spatial  grid.  We  also  need  the  relation  between  the  displacement 

I _ I  ' 

current  and  the  electric  field;  namely,  p  =  e(a>)E  where  e(a>)  is  the  frequency  dependent  permittivity 
and  is  related  to  the  refractive  index  n  by  n2  =  s(co) .  In  the  many  papers  we  have  published  using  the 
FDTD  method  to  solve  scattering  from  irregularly  shaped  particles  (see  references  15-26),  we  have 
worked  only  at  a  single  frequency  since  the  incoming  radiation  was  assumed  to  be  monochromatic. 
Now  when  using  short  pulses,  there  are  a  plethora  of  frequencies  involved  and  so  then  we  must  be  able 
to  use  the  actual  dispersion  relation  for  the  medium  involved.  In  the  case  of  water  in  the  near  IR  (  our 
pulses  are  centered  around  800  nm)  we  will  use  the  actual  measured  absorption  data  of  either  Kuo  or 
Segelstein  in  order  to  apply  the  FDTD  method  to  this  problem,  we  need  to  convert  the  frequency 
dependent  permittivity  to  a  time  dependent  one.  In  particular,  the  temporal  characteristic  of  the  water 
permittivity  should  be  modeled  around  the  overtones  of  the  OH-stretch/scissor  coupled  modes.  These 
are  the  regions  where  the  sharp  changes  in  the  slope  of  the  absorption  coefficient,  denoted  by  a,  occur. 
Another  very  interesting  aspect  of  water  occurs  when  it  is  deuterated  to  form  either  HDO  or  D2O. 
These  forms  will  have  radically  different  absorption  spectra  and  we  can  further  test  our  analyses  by 
studying  these  different  forms  of  water. 

If  we  denote  the  refractive  index  of  water  by  n=nR-inj  and  the  permittivity  by  £  =sR-isIy  where  the 
subscripts  R  and  I  denote  real  and  imaginary  parts  respectively,  then  it  is  easy  to  show  that 
sR  =  n\  -  n]  and  £,  =  2 nRn, .  The  absorption  coefficient  a  is  related  to  ///  by  a  =  4 nnI  /  Avac  where 
A, vac  denotes  the  wavelength  of  the  radiation  in  vacuo. 

In  order  to  get  the  time  dependence  of  water  permittivity,  we  must  Fourier  transfonn  £(co)  .  Actually  it 
is  more  physical  to  Fourier  transfonn  £(a>)  - 1  and  to  obtain  the  susceptibility  kernel  G(  1 )  defined  by  ( 
see  Jackson- ): 


then 


G(t)  = 


2  n 


P  \s{(o)-\\eimd(o 

J  —00 


(3) 


r  r  r  r  r°°  r  r 

D(x,t)  =  E(x,t)  +  G{r)E{x,t  -  r)dv 

J  —00 

r  r  r  r  (4) 

=  E(x,t )  +  y^yG(r)E(x,t  -  t)At 

n  r=° 

Thus  we  see  that  the  displacement  current  at  a  certain  instant  in  time,  depends  on  the  electric  field 
from  earlier  times.  It  should  also  be  noted  that  if  s  is  independent  of  frequency,  then  Eq.  (4)  reduces  to 

I _ I  ' 

p|  =  sE  which  is  what  our  present  FDTD  code  now  uses. 
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Once  we  have  established  the  fact  that  the  FDTD  method  is  a  viable  one  for  solving  this  problem,  then 
we  would  like  to  next  explore  the  possibility  of  using  the  pseudospectral  time  domain  (PSTD)  method 
which  we  feel  will  run  much  faster  than  the  conventional  FDTD  method.  We  have  tested  this  method 
on  some  canonical  scattering  problems  and  found  substantial  reductions  in  computation  time  when 
compared  to  the  conventional  FDTD  method.  With  this  method  we  will  actually  be  able  to  monitor  the 
complete  temporal  evolution  of  the  pulse  shape. 

WORK  COMPLETED 

a)  We  have  demonstrated,  for  the  first  time,  precursor-like  propagation  through  a  bulk  medium 
composed  of  an  organic  dye  solution.  This  is  an  important  first  step  toward  precursor  formation 
and  observation  in  various  other  bulk  media. 

b)  We  simulated  the  propagation  of  a  reconstructed  electric  field  through  the  above  organic  dye 
solution  to  demonstrate  the  formation  of  precursor  features;  theoretical  simulations  and 
experimental  results  were  in  good  agreement. 

c)  In  order  to  characterize  the  optical  pulses  involved  in  precursor  formation,  we  have  developed  a 
novel  and  accurate  pulse  retrieval  algorithm.  This  algorithm  is  a  hybrid  between  the  popular 
graduated  optimization  algorithm  (GOA)  and  the  genetic  algorithm  (GA)  approach;  it  is  able  to 
characterize  sharp  perturbations  in  the  spectral  phase  of  propagating  pulses  more  effectively  than 
other  popular  methods,  and  so  is  ideal  in  describing  the  propagation  of  optical  precursors  through 
optically  dispersive  media. 

d)  We  have  demonstrated,  for  the  first  time,  energy  coupling  between  optical  filament-forming  pulsed 
laser  beams  in  bulk  media.  This  demonstrates  that  filament  formation  and  propagation  can  be 
controlled  after  the  initial  fonnation  of  the  filament  and  has  important  consequences  in  atmospheric 
and  oceanic  remote  sensing,  communications,  and  spectroscopy. 

RESULTS 

We  first  present  the  results  obtained  from  observation  of  precursor  propagation  through  a  bulk  medium 
composed  of  an  organic  dye.  The  dye  itself  (NP800,  Exciton)  was  selected  for  its  sharp  absorption 
resonance  in  the  near  infrared  portion  of  the  electromagnetic  spectrum,  with  the  peak  absorption 
occurring  around  800  nm  with  a  full  width  at  half  maximum  (FWHM)  of  about  100  nm.  The  calculated 
refractive  index  of  the  dye,  obtained  via  the  Kramers-Kronig  relationship,  shows  strong  anomalous 
dispersion  in  the  region  around  800  nm.  This  strongly  peaked  absorption  and  anomalous  dispersion  are 
requirements  for  precursor  formation. 

The  dye  was  dissolved  in  1-1-2-2-tetrachloroethane,  with  a  final  dye  concentration  of  10  g/L.  A  thin 
layer  of  the  solution  was  sandwiched  between  two  microscope  cover  slips,  each  about  190  um  in 
thickness;  surface  tension  held  the  dye  solution  between  the  slips,  and  the  whole  was  placed  in  the  path 
of  the  beam.  The  laser  pulses  were  generated  by  a  femtosecond  laser  system  (Rainbow,  Femtolasers) 
and  had  a  central  wavelength  of  800  nm,  a  temporal  duration  of  9  fs  FWHM,  and  pulse  energy  of  4  nJ. 
The  low  energy  ensures  that  our  experiments  were  completed  within  the  linear  regime  of  the  light- 
matter  interaction. 

We  first  measured  the  autocorrelation  trace  of  the  initial  pulse;  comparison  with  the  autocorrelation 
pulse  after  propagation  through  the  organic  dye  solution  thus  provided  sufficient  information  for  the 
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identification  of  precursor  formation.  Reconstruction  of  the  pulse  electric  field  from  the 
autocorrelation  traces  was  completed  with  the  hybrid  genetic  algorithm  described  below;  the  initial 
electric  field  and  the  electric  field  after  propagation  through  the  organic  dye  solution  are  shown  in 
Figs.  4  and  5,  respectively. 


Fig.  4  The  reconstructed  electric  field  of  the  initial  femtosecond  laser  pulse  as  provided 
by  the  hybrid  genetic  algorithm.  Reference  peaks  A,  B,  and  C  are  labeled. 


1.0 
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Fig.  5  The  reconstructed  electric  field  of  the  femtosecond  laser  pulse  after  propagation 
through  the  organic  dye  solution  as  provided  by  the  hybrid  genetic  algorithm.  Reference 
peaks  A,  B,  and  C  are  labeled.  Notice  that  magnitude  of  the  peaks  have  attenuated  at 
differing  rates  as  a  result  of  propagation  through  the  dye  solution. 


The  simulated  autocorrelation  traces  using  the  reconstructed  electric  fields  were  in  good  agreement 
with  the  measured  traces.  Theoretical  numercial  simulations  propagating  the  intitial  electric  field  were 
then  completed  via  the  generalized  pulse  propagation  equation  to  visualize  the  shape  of  pulse  features 
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as  a  function  of  propagation  distance.  The  peak  magnitude  of  the  pulse  features  labeled  as  A,  B,  and  C 
in  Figs  4  and  5  were  tracked  over  propagation  distance  x  to  the  thickness  of  the  organic  dye  solution  d. 
The  results  are  shown  in  Fig.  6.  One  can  see  the  differential  attenuation  of  these  pulse  features  as  a 
function  of  propagation  distance.  Such  behavior  is  characteristic  of  precursors. 


1.2 
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Fig.  6  Simulated  attenuation  behavior  of pulse  features  A,  B,  and  C  from  Figs.  1  and  2 
over  the  thickness  of  the  organic  dye  solution.  Notice  the  varying  attenuation  coefficients; 
such  behavior  is  characteristic  of  precursor  propagation. 


Understanding  of  precursor  formation  and  propagation  in  bulk  media  is  important  in  realizing 
sophisticated  methods  of  utilizing  pulsed  laser  beams  for  communication  and  remote  sensing.  In 
addition,  precursors  themselves  may  perhaps  find  applications  in  characterization  of  material 
properties. 

We  now  briefly  present  the  hybrid  genetic  algorithm  developed  to  retrieve  the  above  electric  fields 
from  the  autocorrelation  traces. 

The  graduate  optimization  algorithm  (GOA)  was  able  to  scan  the  chirp  parameters  of  the  spectral 
phase  of  ultrashort  pulses.  Combined  with  the  measured  spectral  intensity,  spectral  phase  can  be 
coarsely  retrieved  by  matching  retrieved  and  measured  second  harmonic  interferometric 
autocorrelation  (SHIAC)  within  a  calculation  time  of  seven  seconds. 

The  genetic  algorithm  (GA)  was  then  incorporated  into  the  retrieval  and  employed  to  make  fine 
adjustments  of  both  individual  spectral  phases  and  the  chirped  parameters  to  further  reduce  the 
difference  between  retrieved  and  measured  SHIAC.  The  calculation  took  several  minutes  or  more  than 
ten  minutes  depending  on  how  many  iterations  were  used  in  the  calculation. 
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Fig.  7  shows  the  GOA+GA  combined  algorithm  retrieved  a  spectral  phase  that  matches  with  measured 
SFIIAC  perfectly,  (a)  shows  the  retrieved  SFIIAC  from  GOA  only,  (b)  shows  the  retrieved  SHIAC 
from  GOA+GA.  There  is  a  perfect  match  between  retrieved  and  measured  SHIAC.  (c)  shows  the 
retrieved  phase  from  GOA  (blue  curve)  and  GOA+GA  (red  curve),  (d)  shows  the  evolution  of  rms 
errors  between  retrieved  and  measured  SHIAC. 


0  50  100  150  200 
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Fig  7.  (a)  Retrieved  SFIIAC  from  GOA  only,  (b)  Retrieved  SHIAC  from  GOA+GA.  (c)  Retrieved 
phase  from  GOA  (blue  curve)  and  GOA+GA  (red  curve),  (d)  Evolution  of  rms  errors  between 

retrieved  and  measured  SHIAC. 

Fourier  transform  of  the  above  spectrum  with  corresponding  spectral  phase  yields  the  pulse  shown  in 
Fig.  8. 


Fig.  8  Retrieved  electric  field  from  the  hybrid  genetic  algorithm  shown  in  Fig.  7. 
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GOA+GA  also  works  well  with  more  complex  pulse  shapes,  as  shown  in  Fig.  9.  (a)  shows  the 
retrieved  SHIAC  from  GOA  only,  (b)  shows  the  retrieved  SFIIAC  from  GOA+GA.  There  is  a  very 
good  match  between  retrieved  and  measured  SFIIAC.  (c)  shows  the  retrieved  phase  from  GOA  (blue 
curve)  and  GOA+GA  (red  curve),  (d)  shows  the  evolution  of  rms  errors  between  retrieved  and 
measured  SFIIAC. 
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Fig.  9  (a)  Retrieved  SHIAC  from  GOA  only,  (b)  Retrieved  SHIAC  from  GOA+GA.  (c)  Retrieved 
phase  from  GOA  (blue  curve)  and  GOA+GA  (red  curve),  (d)  Evolution  of  rms  errors  between 

retrieved  and  measured  SHIAC. 

Fourier  transfonn  of  the  spectrum  and  corresponding  spectral  phase  yields  the  pulse  shown  in  Fig.  10. 
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Fig.  10  Retrieved  electric  field  from  the  hybrid  genetic  algorithm  shown  in  Fig.  6. 
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Hence  we  have  shown  that  the  developed  genetic  algorithm  is  able  to  accurately  reconstruct  pulses 
which  possess  complex  and  sharply  varying  spectral  phases,  and  thus  enjoys  an  advantage  over  other 
less  intensive  approaches.  Reconstruction  of  complex  pulse  shapes  is  important  in  understanding  their 
propagation  through  media  which  possess  highly  dispersive  and/or  absorptive  refractive  indices,  and 
the  developed  algorithm  is  therefore  a  valuable  tool  in  studying  their  pulse  propagation  behavior. 

Having  reported  investigations  into  linear  phenomena,  we  now  present  experimental  results  involving 
nonlinear  phenomena.  Optical  filaments  are  complex  light  pulses  resulting  from  intensity-dependent 
Kerr  lensing  and  self-focusing;  the  Kerr  focusing  is  balanced  by  plasma  defocusing  effects,  thus 
allowing  an  ultrashort  laser  pulse  to  propagate  large  distances  without  diverging.  They  have  promising 
applications  in  sophisticated  techniques  for  atmospheric  and  oceanic  remote  sensing,  imaging, 
communication,  and  spectroscopy. 

One  important  parameter  in  these  applications  is  control  of  the  propagating  filament.  This  may  be 
achieved  through  beam-coupling  between  two  filament-forming  pulsed  beams.  Such  coupling  has 
already  been  demonstrated  in  air  (A.C.  Bernstein  et  al.  (2009),  Y.  Liu  et  al.  (2010));  we  have 
demonstrated,  however,  the  first  instance  of  energy  coupling  between  filament-forming  beams  in  bulk 
media.  Moreover,  we  report  an  undocumented  phenomenon  in  which  the  energy  transfer  oscillates  on 
an  increment  in  relative  pulse  delay  equal  to  the  time  period  of  the  incident  radiation. 

We  used  a  Ti:sapphire  based  laser  system  consisting  of  an  oscillator  (Mira,  Coherent)  and  an  amplifier 
(Legend,  Coherent:  800  mn  center  wavelength,  35  fs  transform-limited  pulse  duration,  1  kHz  repetition 
rate,  1  mJ  pulse  energy).  Two  identical  chirped  pulses  were  generated  and  crossed  at  a  small  angle 
within  a  glass  cell  containing  methanol.  We  chose  methanol  as  the  bulk  medium  because  of  its 
transparency  at  the  wavelengths  used;  it  is  also  an  important  first  step  toward  realizing  energy  coupling 
in  other  media.  The  energy  of  the  pulses  was  about  2  uJ,  enough  to  form  filaments  after  passing 
through  a  lens  with  a  focal  length  of  5  cm.  The  pulse  energy  was  then  varied  as  an  experimental 
parameter.  The  resulting  filament  energies  were  then  measured  with  a  spectrometer  (USB2000,  Ocean 
Optics).  As  we  used  chirped  pulses,  the  relative  delay  between  the  pulses  determined  the  direction  and 
amount  of  energy  transfer;  this  was  controlled  by  passing  one  beam  through  a  translation  stage 
(ESP301,  Newport).  We  refer  to  the  two  beams  as  the  fixed  and  translation  stage  beams,  respectively. 

Results  are  shown  in  Fig.  11.  Six  curves  are  shown  for  various  pulse  energies.  It  is  important  to  notice 
that  lower  pulse  energies  show  energy  transfer  consistent  with  A.C.  Bernstein  et  al.  (2009),  while 
higher  pulse  energies  show  energy  transfer  consistent  with  Y.  Liu  et  al.  (2010).  This  is  interesting  in 
that  behaviors  documented  for  filament  energy  transfer  in  air  are  also  present  in  bulk  media. 
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Fig  11.  Energy  transfer  curves  for  various  beam  powers  in  methanol.  Lower  powers  show  energy 
transfer  consistent  with  A.C.  Bernstein  et  al.  (2009),  while  higher  power  levels  show  energy  transfer 

consistent  with  Y.  Liu  et  al.  (2010). 

The  curves  shown  in  Fig.  8  appear  “noisy”;  however,  it  can  be  shown  that  the  above  curves  actually 
oscillate  at  the  same  frequency  as  the  incident  radiation.  This  is  shown  in  Fig.  12. 
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Fig.  12  (a)  Energy  transfer  curve  for  filament-forming  pulses  propagating  in  methanol,  (b)  Zoomed 
in  region  of  the  same  curve  in  (a)  with  superimposed  sinusoid  having  time  period  of  2.6  fs,  the  same 

as  the  incident  radiation. 
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(a)  shows  an  energy  transfer  curve  for  two  filament  forming  beams  in  methanol,  while  (b)  shows  a 
zoomed-in  region  of  the  same  curve  with  a  superimposed  sinusoid  have  the  same  time  period  as  the 
incident  radiation. 

These  results  indicate  that  very  fine  control  of  filaments  may  be  achieved  through  two-beam  coupling. 
Many  more  applications  for  these  results  may  soon  be  realized. 

IMPACT/APPLICATION 

The  work  we  are  doing  on  nonlinear  filament  formation  may  lead  to  a  very  new  and  novel  way  to  do 
hyperspectral  remote  sensing  and  underwater  communication. 

TRANSITIONS 

If  the  results  of  this  research  uncover  novel  ways  to  propagate  acoustic  signals,  then  we  will  have  a 
scheme  to  measure  the  speed  of  sound  in  water  very  quickly  and  accurately  from  helicopters  or  aircraft 
covering  large  regions  of  the  ocean. 

RELATED  PROJECTS 

We  use  the  results  from  our  research  on  the  RaDyO  project  to  give  us  the  IOP’s  for  water  we  use  in 
our  analysis  for  short  pulses. 
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